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a. 

Abstract. We analyze in detail the Solid-On-Solid model (SOS) for growth processes 
on a square substrate in 2 + 1 dimensions. By using the Markovian surface properties, 
. we introduce an alternative approach for determining the roughness exponent of a 

special type of SOS model-the Restricted-Solid-On-Solid model (RSOS)- in 2 + 1 
dimensions. This model is the SOS model with the additional restriction that the 
height difference must be S = 1. Our numerical results show that the behaviour of 
the SOS model in 2 + 1 dimensions for approximately S > S x ~ 8 belongs to the 
two different universality classes: during the initial time stage, t < t x it belongs to 
O \ the Random-Deposition (RD) class, while for t x < t <C t sat it belongs to the Kardar- 

Parisi-Zhang (KPZ) universality class. The crossover time (i x ) is related to S via a 
power law with exponent, r\ — 1.99 ± 0.02 at la confidence level which is the same 
as that for 1 + 1 dimensions reported in Ref. |S]. Using the structure function, we 
compute the roughness exponent. In contrast to the growth exponent, the roughness 
, exponent does not show crossover for different values of S. The scaling exponents 

of the structure function for fixed values of separation distance versus S in one and 
two space dimensions are £ = 0.92 ± 0.05 and £ = 0.86 ± 0.05 at la confidence level, 
| respectively. 
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1. Introduction 

Surface growth processes, especially the formation of thin film deposits, have 
been studied using various approaches in complex systems and statistical analysis 
[II El El HI EJ EJ E] . The factors which control surface growth phenomena have immense 
phase space. Consequently, to be able to analyze these phenomena one needs to make 
many assumptions, which can lead to results that are unreliable. Combining insights 
from computational simulation and simplified analysis will likely give better results. 
It is well-known that the understanding of phenomena such as advances of bacterial 
colonies, electrochemical deposition, flameless fire fronts and molecular-beam-epitaxial 
growth is of considerable importance in the control of many interesting growth processes 
in industries [U El E]. The simplest surface growth model is the so called statistical 
deposition model [TJ [TO]. Some models proposed to explore growth surfaces, such as 
the Family model [11] , Ballistic Deposition(BD) model [T2"l [13] and the Eden model |14j . 
are able to account for many of the properties of some real systems. For example, the 
BD and Eden models can accurately simulate vapor deposition and biological growth. 
However, these models tend to ignore the microscopic details of the interfaces , and 
cannot provide accurate scaling exponents. In addition many fractal features of real 
systems remain unexplained [9J [THJ [T6] . To solve these problems, one should modify the 
above models. 

The Solid-On-Solid model (SOS) is more suitable to describe a real surface's 
properties than those models described above PITH [El EE]- This growth model does not 
exhibit strong corrections to scaling and consequently allows us to determine accurate 
values of scaling exponents 0[T71[IE]. The Restricted-Solid-On-Solid (RSOS) model( a 
modified version of the SOS model), proposed by Kim et al. |17j , is most important due 
to its wide applicability, such as for surface roughening modeling via exothermic catalytic 
reactions on the substrate [8]. Various aspects of the Solid-On-Solid model for surface 
growth have been studied: the effect of long-range elastic interactions [20], growth 
processes with correlated noise [21] , phase transitions as a function of temperature-like 
parameters |22j, the (OOl)-surface morphology of GaAs annealed at fixed temperature 
and pressure, the well explained by annealed version of the RSOS model [231 121] ■ 
Crossover from random to correlated regime [251 [26] , relaxation to steady states [27], 
distribution of local configurations for finite values of S [8], Markov analysis [28], the 
effect of hopping in various local growth rules on the linear and nonlinear fourth-order 
dynamical growth equation [29], growth model in higher dimensions |30j and, more 
recently, the growth on fractal substrates based on the SOS model [51] . has also been 
addressed in the literature. 

As mentioned in many previous studies, it is believed that the RSOS model belongs 
to the Kardar-Parisi-Zhang (KPZ) universality class in the continuum limit [32J |33] . The 
KPZ equation is one of the most important phenomenological theories in which time 
evolution of the interface has been characterized by the height function h(f, t) at position 
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f and time t. The governing equation is given by 

= vV 2 h(r,t) + ^[Vh(r,t)] 2 + K(f,t). (1) 

Here v and A represent the surface tension and the excess velocity respectively, while 
K(f, t) is a Gaussian noise with zero mean and co- variance 
(K(r, t)K(r', t')) = D5 d (r- ?)5(t - if), 

where d is the dimension of the substrate, and D is the noise intensity [H [12]. The 
interface width reads as: 

1/2 v 



W(L d ,t) 



1 2>tf*) - M*)] 2 



L d 



(2) 



This characterizes the roughness of the interface, for growth in a substrate of length L, 
and h(t) is the spatial average of height at time t. For short times, the interface scales 
as follows: 

W(L,t)fatP (3) 

where /3 is called the growth exponent. For long times, a steady state is attained and 
the width is saturated as follows: 

W sat (L,t)^L a . (4) 

Here a is the roughness exponent. Equations ([3]) and (j3J) correspond to limits of the 
dynamical relation of the Family and Vicsek ansatz: 

W{L,t)*L«f(^. (5) 

The dynamical exponent, z — % , characterizes the crossover from the growth regime to 
the steady state. The exact scaling exponents are known in d = 1, but no exact value 
has been obtained in two or more dimensions [IE] . Many discrete models fall into the 
KPZ class, such as the RSOS model [TTl [18] and Ballistic Deposition BD [12]. Most of 
the reported values of a range from a = 0.37 to a = 0.40 [TT], [TBI G2E ESI EZ], confirmed 
by numerical solutions of the KPZ equation [381 E3 HQ]. 

The competition between different growth mechanisms during particle deposition, 
as well as phase transitions which are very often observed in many real growth processes, 
has been investigated in many studies [25], [26] . Recently, it has been confirmed that there 
exists a crossover between the Random Deposition and KPZ classes at the initial growth 
stages for all values of the height restriction parameter between nearest neighbours for 
the SOS model in 1 + 1-dimension [8]. Here we are interested in investigating the 
possibility of the existence of crossover in the SOS model in 2 + 1 dimensions. In 
addition, we give a new approach to determine the roughness exponent using Markovian 
properties of surfaces. 

The rest of this paper is organized as follows: in section 2, we introduce the 
Markovian surface, and by using the characteristic function, the roughness exponent is 
calculated. The SOS model for finite values of S is numerically investigated in section 3. 
Crossover in the growth mechanism and corresponding properties are also investigated 
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Fig.l A typical trajectory from point r = (0,0) to r = (n x ,n y ). 

2. Markovian Surface 

The Markovian surface is one of several models to represent multi-level (stepped) 
crystalline surfaces. In this model, it is assumed that the steps have only a mono-atomic 
height. Displacement through any steps may be upward or downward, each occuring 
with equal probability. Let 7 be the probability of meeting an atom displaced vertically 
either upward or downward in going from any lattice site to an adjacent one. That is, 
the probability of encountering a step (Ah = ±1) while the corresponding probability 
for a lateral walk, namely Ah = (h is the height of the surface), is equal to 1 — 7. 
Since every displacement or step occurs independent of any other, the step surface is 
mapped to the path of a Markovian chain |41j . For the Markovian chain or random walk 
model, there exist three choices for the displacement at each walk: an upward walk with 
a probability 7/2, a downward walk with a probability 7/2 and a lateral walk with a 
probability 1 — 7 [H]. As mentioned in the introduction, here we rely on the Markovian 
surface to explore the scaling exponent of the RSOS growth model. To this end, we 
introduce the characteristic function defined as the Fourier transform of the probability 
distribution function, P(Ah(r)), with respect to Ah(r) = h(r) — h(0) after saturation 
time, as 

Z d (X,r} = (e lXc[h ^- hm ) (6) 

where c is the unit of step variations, which is equal to one in the Markovian surface 
and RSOS model. The height difference, h(r) — h(0) can be represented as the sum of 
the height differences between successive sites from r = to r = na in one-dimension 
and to r = J(n x 2 + n y 2 )a in two- dimensions (a is lattice unit). In one-dimension we 
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have 

h(na) - h(0) = ^ [h(ia) - h((i - I) a)] (7) 

8=1 

For the RSOS model in 2 + 1 dimensions, the height difference between any sites with 
coordinates (n x ,n y ) and its nearest neighbour sites with coordinates (n x ±l,n y ) and 
(n x ,n y ± 1) is ±1. To calculate the characteristic function, we should move from point 
r = (0,0) to f = (n x ,n y ) in different paths like path A as shown in Figure (1). So the 
vector sum of the trajectories within path A gives the vector r. Due to isotropy and 
homogeneity of the surface, the probability 7 is similar for each step. Consequently, the 
Z d {X,r) can be written as follows: 

Z d (\,r)\ A = (e iXc M f >-m))\ A 

_ ^J,Xc[h(n x ,n y )—h{n x —l,n y )\\ 

x fJ,Xc[h{n x — l,n y — 1)— h(n x — 2,n y )]\ 
x ^i\c[h(n x -2,n y )-h(n x -2,n y -l)}^ 

(8) 

The number of different paths from point r = (0, 0) to r = (n x , n y ) is (n x + n y — 2) with 
the condition n x , n y > 1 and the number of paths from point r = (0, 0) to r = (n x , n y ) 
is N = Cl x x n + ny -i = {n x n y )\/{n x n y - n x - n y + l)\{n x + n y - 1)!. 
Finally the characteristic function Zd(X,r) can be written as 

Z d (\,r) = Nf n * +n y- 2 . (9) 

Here 

f — / e i\c[h(n x ,n y )—h(n x —l,n v )]\ 

dAh e lXcAh(ff) P{Ah{f)) 
= 1 - 7[1 - cos(Ac)] (10) 
Therefore, the characteristic function for the RSOS model is 

Z d (X,r) = N{l- 1 [l- C os(\c)}} {n * +n «- 2) (11) 
For r — >■ regime and small 7, equation (fTlj) becomes: 

Z d (X,r) = e 7[l-co S (Ac)](n a; +n !/ -2) ^ 

On the other hand, the above equation can be expressed [16] as 

Z d (\,f) = e-^ 2a (13) 
By comparing equations ([12]) and ( IT3|) . one finds: 

- 2Q = (*4+«5)° 

~K + n y -2) (14) 

In order to determine the roughness exponent, a, we have to compute the following 
scaling relation for small f: 

n x + n y - 2 ~ (n x 2 + n y 2 ) a (15) 
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By using numerical calculations and averaging over small values of r, the roughness 
exponent can be read as 0.39 ± 0.03 in la confidence level, which is in good agreement 
with the previous results [T7| [T8| [35j [36| [37]- 




J ! ! 1 ! ! ! 1 ! ! ! 1 

10' l(f 10* 

t 



Fig. 2 Log-log plot of interface width versus time for various values of parameter S. Here the 
system size is L x L = 4096, and ensemble averaging has been done over 500 independent 
runs. Solid lines show the scaling behavior of RD and KPZ models in 2 + 1 dimensions with 
slopes (3 = 0.5 and /3 = 0.25, respectively. 

3. Simulation of Solid on Solid growth model in 2 + 1 dimensions 

In the previous section, we dealt with the Markovian surface and calculated 
roughness exponent for restricted solid-on-solid model in a new way. In this section 
we simulate this surface growth model in 2 + 1 dimensions for finite values of parameter 
S on the square lattice with length L. The growth process of the SOS model can be 
described by the following steps: 

Step 1: Select a site randomly e.g. site f : (n x ,n y ) among all L 2 sites. 

Step 2: Then all of the following conditions should be satisfied to increase the height of 

mentioned site: 

I) | [h(n x , n y ; t) + l]- [h(n x - 1, n y ; t)] \ < S 

II) | [h(n x , n y ; t) + 1} - [h(n x + 1, n y ; t)} \ < S 

III) | [h(n x , n y ; t) + l]- [h(n x , n y - 1; t)] | < S 

IV) | [h(n x , n y ; t) + 1] - [h(n x , n y + 1; t)} \ < S 

V) Otherwise, do nothing. 
Step 3: Repeat the above tasks. 

To reduce the errors due to the substrate's boundaries, we use periodic boundary 
conditions during particle deposition. Each time step is defined as the number of 
particles needed to fill the surface average , which is equal to L x L [Tfl HE]. The 
log-log plots of the W(L 2 ,t) versus time scale for SOS model with different values of S 
are shown in Figure (2). The slope of these diagrams, for the initial time scale of the 
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growth process, gives the growth exponent, /3, of the model. Figure (2) demonstrates 
that there exists crossover for t -C t sat in the log-log plot of interface width versus growth 
time. To investigate the crossover behavior of interface width at small time for various 
values of parameter S, we define a fluctuation function as: 

MS) 

MS)= E \W(S,t)-W$ he (S,t)\ (16) 
t 

where W^ hc (S, t) ~ and the symbol o stands for RD and KPZ and corresponding 
exponents are 0.50 and 0.25, respectively. It is worth noting that RD and KPZ are 
the most relevant universality classes for SOS growth model reported in [8j [T7J [18] . As 
Figure (3) shows, by increasing the value of the restricted parameter S, deviation(5S') 
from KPZ(RD) class at the initial growth stage increases(decreases). After enough time, 
the interface width of SOS model in 2+1 dimensions will be saturated. The value of S for 
the intermediate value of fluctuation function (equation ( fl6|) ) is S x = 8. Subsequently 
one concludes the SOS model in 2 + 1 dimensions for approximately S > S x and before 
saturation epoch belongs to the two different universal classes: at the very early growth 
stage, it belongs to the RD class and at intermediate time scales, t x < t <C t sat , it tends 
to the KPZ universality class and is affected by the restrictions on the height differences 
or behaves like the RSOS model with S = 1. On the other hand, for S < S x , the SOS 
model only belongs to one class, namely, KPZ class, before the saturation of its interface 
width. The mathematical form of this dynamic for S > S x can be read as 

W(S,t) = W(S)f(—) (17) 



in which 



t 



,0.50 

(-) ~< )7(0.25 (18) 



where r x is the crossover time scale indicating transition from random deposition 
universality class to RSOS class which in principle depends on S. The quantity W(S) 
just depends on S. Recently Ching-Chun Chien et al.[8] have shown that the growth 
exponent, 0, for SOS model in 1 + 1 dimensions during the very early stage of surface 
growth is independent of S. Their results confirmed that SOS model in 1 + 1 dimensions 
for t < t x belongs to RD class and then crosses over to the KPZ universality class for 
t x < t <C t sat , while the extension of their results to 2 + 1 dimensions demonstrates 
that for S < S x , it belongs to only one universality class for t t sat (Figure (3)). 

In order to interpret the crossover in the interface width function (see Figure (2)), 
we refer to the correlation function, C(r,r') defined by 

C(r,P) = (\h(r,t) - h][h(?,t) - h}) (19) 

For an isotropic surface, we can define the normalized correlation function as follows: 

C(\* ?l) _ (ihir + tt) -h}{h(r,t) -h}) 
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The correlation functions for various values of S at the early time scale, t <C t sat , are 
shown in Figure (4). The correlation length scale, over which the correlation function 
reaches 1/e of its maximum, decreases as S increases. Moreover, as S increases the effect 
of randomness in the particle deposition process decreases. This can be explained as 
follows : for larger values of S, during the very early stage of surface growth, the height 
of a typical site is not affected by its neighbouring sites due to the restriction constraint 
embedded in the rule of its deposition. Consequently, one expects the correlation length 
of height to decrease as S increases at the very early stage of growth. Therefore, the 
SOS model with infinite S reduces to the RD model during the very early growth stage. 
Figure (5) indicates the log-log plot of time in which the SOS model crosses over from 
RD class to KPZ class, r x versus S. This confirms the scaling behaviour of transition 
time versus height difference. The slope of this plot is rj = 1.99 ± 0.02 at la confidence 
interval. This value is in agreement with results in 1 + 1 dimensions given by Chih-Chun 
Chien et al. [8]. It may be stated that exponent has the same values in 1 and 2-space 
dimensions, while the growth and roughness (see below) exponents depend on space 
dimension. 




Fig.3 The function of A(S) for /3 KPZ = 0.25 (filled circle symbol) and (3 RD = 0.50 (filled 
square symbol) at t <C t sat . 
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Fig. 4 Correlation functions for SOS model for different values of height restriction 
parameters. 



i — i — i i i i 1 1 1 1 




Fig. 5 The crossover time scale as a function of parameter S. Solid line indicates the scaling 
function with exponent, n = 1.99. 
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Fig. 6 Log-log plot of height correlation function versus separation distance for various values 
of S. Solid line corresponds to scaling function with exponent, a = 0.40. 

Another interesting parameter is the roughness exponent, a, which is also given via 
the structure function for each S as follows: let the structure function be defined as 



G(r,r';t) = ([h(f,t) - h(r' } t)Y) 



2\V2 



(21) 



In the homogeneous and isotropic case G depends on time scale, t and separation 
distance, I — \r — r'\. In addition, in principle this function may depend on S, so 
in the most general case, the structure function reads as: 



G(l,S;t) 



(22) 



a//3. The roughness saturates after a sufficiently long time; consequently 
behaves as: 



where z - 
G(l,t>t sat ) 

G(l, Sfixed, t > t sat ) = G(l, Sfixed) ~ 1° (23) 

Therefore, the slope of log-log plot of G(l, 5fi xe d) versus I for small separation distance in 
the saturation regime gives the roughness exponent, a. Figure (6) shows the structure 
function for different values of S in the saturation regime. The slopes of G(l, Sfixcd) for 
small I are all equal to 0.40 ± 0.01 at la confidence level which is the same as that 
determined for the KPZ model in 2 + 1 dimensions. This confirms that the SOS model 
belongs to the KPZ universality class during the late growth stage. Figure (7) shows the 
log-log plot of G(l, S) for small separation distance and fixed / versus S. It demonstrates 
a scaling behavior for structure function in the saturation regime for small and fixed / 
as a function of S. Its scaling exponent is equal to £ = 0.86 ± 0.05, at 68% confidence 
level. We introduce a new scaling function which gives the relation between I and S 
after saturation epoch: 



G{l,S) 



(24) 
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Table 1. Values of the scaling exponents of SOS growth model in 1 + 1 and 2 + 1 
dimensions. 



Dimension 









V 


a 




1 + 1 


All S 

t < t x 
t > t x 


[8] 
0.5 
0.33 




2.06 [8] 


0.5 


0.92 + 0.05 


2 + 1 


s <s x 
s> s x 
s>s x 


0.25 

t < t x 
t > t x 


0.5 
0.25 


1.99 + 0.02 


0.40 + 0.01 


0.86 + 0.05 



where z' = j which is a new dynamical exponent. We also examined the height-height 
correlation function in 1 + 1 dimensions for various values of S. Our results confirm 
a scaling behavior with the exponent equal to £ = 0.92 ± 0.05. Table (PQ) reports all 
the most relevant exponents determined in this paper as well as those given in ref. [8]. 




i i i i i i 

1 2 3 4 5 6 7 S 

S 



Fig. 7 Log- log plot of G{1 = fixed) versus S in 2 + 1 dimensions. Solid line shows the scaling 
function with slope £ = 0.86. 

4. Summary and conclusion 

In this paper, we explored some scaling properties of the Solid On Solid model, one 
of the well-established algorithms for surface growth in 2 + 1 dimensions. The scaling 
properties of the model give deep insight into its universality and statistical classification. 
To this end, we have proposed a new method for finding the roughness exponent for 
one of the special classes of the SOS model: the so-called RSOS model. In the RSOS 
model, the value of near-neighbours height difference is restricted to S = 1. This new 
method uses Markovian surfaces properties. By introducing a characteristic function 
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we computed the scaling behavior of the site's numbers with respect to the various 
values of the displacement vector. The value of the roughness exponent derived by the 
Markovian approach is a = 0.39 ± 0.03 at 68.3% confidence interval. Our result for the 
roughness exponent is in good agreement with those obtained from direct simulation 
of the RSOS model performed in earlier works pH HHJ [35J |36j [37]. We simulated the 
growth surface according to the SOS algorithm explained in section 3 for finite values of 
S. The interface width as a function of time has a crossover time scale for approximately 
S > 8 which confirms two universality classes for the SOS model at the very early stage 
of its growth. The SOS model for S > 8 at t < r x falls under the random deposition 
universality class and thereafter tends to the KPZ universality class for t x < t <C t sat . 
To make these crosses over more obvious, we investigated the normalized correlation 
function. Our results indicate that by increasing the height difference parameter at 
t <C t sat , the correlation length scale is decreased (see Figure (4)). This is clearly due 
to the restriction constraints which eliminate the effect of neighbors on the memory 
of particle deposition. With increasing growth time, particle depositions are affected 
by the restriction rule; consequently, the SOS model with infinite S tends to the KPZ 
growth model before the saturation time stage. According to Figure (5), the scaling 
exponent of r x versus S is rj = 1.99 ± 0.02 which is in agreement with that obtained in 
1 + 1 dimensions [8], while the growth and roughness exponents depend on the space 
dimension. In order to compute the roughness exponent, we used a structure function. 
The slope of this structure function versus distance separation, I in the log-log scale for 
various values of parameter S is a = 0.40 ± 0.01 at 68% confidence level. This quantity, 
in contrast to the growth exponent, did not show crossover behaviour in the log- log plot 
of structure function as a function of separation distance (see Figure (6)). The value 
of the roughness exponent also confirms our result regarding the universality class of 
the SOS model at longer times. Based on Figure (7), we found that the height-height 
correlation function, G(l) for small and fixed I versus S indicates a scaling behavior with 
exponent £ which is equal to 0.86 ± 0.05 and £ = 0.92 ± 0.05 at la level of confidence 
for 2 + 1 and 1 + 1 dimensions, respectively. Using the roughness and £ exponents we 
introduced a new dynamical exponent as z' = at/£. 
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